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THE EQUATIONS OF MOTION IN GENERAL RELATIVITY THEORY 
AND THE ACTION PRINCIPLE 


By L. InFELD’ 
(Received December 17, 1956) 


This paper deals with the formulation of the equations of motion for p bodies, assum- 
ing: 1. the field equations of General Relativity Theory, 2. the fact that the energy- 
-momentum tensor is different from zero only along certain world lines; this means that 
this tensor depends linearly on Dirac‘s 6-functions. 

It can be shown that from these simple assumptions the proper Lagrangian can be 
found almost without tedious calculation. It will give the equations of motion up to the 
post-Newtonian order. 

Also the connection between these equations and those of the “geodetic lines“ is 
elucidated. 


§ 1 Historical introduction 


The history of the problem of motion in general relativity theory usually, and 
rightly so, starts with the paper written by Einstein and Grommer in the year 1927. 
This paper showed for the first time that the equations of motion for a testing particle 
(with mass m —'0) need not be added to the field equations — that they can be de- 
duced from the relativistic field equations. For many years afterward, Einstein, and 
then Einstein with his collaborators, tackled the problem of the motion of two particles; 
to find whether their equations of motion can be deduced from the field equations. 
The answer to this question was given in the paper by Einstein, Infeld and Hoffmann 
in 1938 in which the two-body problem was solved for the first time. 

Independently of us and a little later Fock (1939) also deduced the equations of 
motion, though only the Newtonian ones, from the field equations. Later Papapetrou 
(1951) simplified his procedure and deduced the post-Newtonian equations of motion, 
after Petrova (1949) had done the same thing on the basis of Fock's theory. Petrova's 
and Papapetrou's results were the same as ours. ; 

What are the essential similarities and differences between Einstein’s theory, 
especially as formulated in the two later papers (Einstein-Infeld 1939 and 1949) and 
the Fock-Papapetrou papers? | 

Without going into detail we can say in general that they have one essential idea 
in common but two different ideas, of which I regard only one as essential. 
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The idea common to both sets of papers is the approximation method about which I 
shall say more later. 

Now, about the first difference: 

In the school represented by Einstein we take as our field equations those of empty 
space, that ist in the usual notation: 


1 
Gop = Rap— PEL RISO; (1.1) 
Einstein always thought that to use 
l 
Gag = Rag — 9 Zap R = (8) 7 Tap (1.2) 


(Тав being the energy-momentum tensor) instead of (1.1) is somehow in bad taste; 
because we do not know in (1.2) what T,, is and we mix a geometrical tensor on the 
left side with a physical tensor on the right side. This was the reason for Einstein's 
long search for a unified field theory in which such a mixture of physics and geometry 
would not appear. 

We know that there is no solution of (1.1) representing spread out continuous matter. 
Therefore by assuming (1.1) we represent matter by means of singularities. The method 
first used by us, consisted in forming certain two-dimmensional surface integrals 
over surfaces enclosing these singularities. The field equations prescribed the laws 
by which the surfaces enclosing the singularities and hence these singularities moved. 
Therefore these laws were — until the post-Newtonian approximation —- deduced 
from the field equations. 

However Fock and Papapetrou consider (1.2) and use definite expressions for Tg- 

This difference does not seem to me to be an essential one for the following reasons: 
we do not know the real distribution of matter. Neither of these methods depicts reality 
properly. The use of our method based on (1.1) means: if the two bodies are a great 
distance apart so that we may assume, approximately, central symmetry of the field 
near one body, then the exact knowledge of density distribution inside the enclosing 
surface is not essential. Outside the enclosing surface, (1.1) is valid. 

I would like to characterise this difference between the two schools by an example 
of a simple situation: that is classical gravitational theory. There we have two kinds 
of equations depending on whether matter is represented by singularities or by a con- 
tinuous distribution. In the first case we have Laplace’s equation 


Ag = 0, (L3). 


---——————————— —J— ———— Quar БЕЕ 


1 Greek indices run from 0 to 3. Latin from 1 to 3. Repetition implies summation. The quadratic 
form for a geodetic coordinate-system is 


ds? — Nap аха dx; Noo = 1; Nom = 0; тп = —ô 
The velocity of light с = 1. 


тт’ 
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in the second case Poisson’s equation 


Дф = Ang. (1.4) 


Now it is more common to write Laplace’s equation, when thinking about its 
spherically symmetric in the form 


where m is the mass and ô is the three-dimensional Dirac д — function. 

Let us try to use the Ó-functions consistently in general relativity theory. Here, 
too, we have, in the first approximation, using (1.1) for empty space, the Laplace 
equation which we solve by a central symmetric solution. This means that a solution 
of (1.1), in the first approximation, valid everywhere, corresponds to that of (1.2) 
where the energy momentum tensor Тв is proportional to Dirac's 6-function. Thus 
the use of (1.2) with Тв proportional to дз, corresponds exactly to our previous 
considerations of (1.1) with singular solutions; just as equation (1.5) is only a different 
form of (1.3). 

Indeed the use of (1.2) with Тв proportional to дз) -functions tremendously 
simplifies the entire deduction of the equations of motion. This simplification was 
achieved in my paper of 1954, but the entire procedure can still be clarified and 
simplified. In fact I shall present here the derivation of the post-Newtonian equations 
of motion with almost no tedious calculations?. 

Summarizing, I may say, that though I believe in Einstein's idea of not using 
the energy momentum tensor, yet I was unfaithful to it, because I used the energy- 
momentum tensor as proportional to the ó(,)-function, to express the singular solutions. 

The next difference between Fock's and Einstein's school is more essential. 

Fock, Petrova and Papapetrou use the harmonic coordinate system, that is the 


four equations: 


1 
ah __0\2 ров 

b а) 3 6%], = с (1.6) 

As а matter of fact, Fock considers the choice of this coordinate system to be 
extremely important, claiming that its addition to the gravitational equations (plus 
some conditions at infinity) restricts the coordinate system up to.a Lorentz trans- 
formation. Thus, for Fock, the choice of the harmonic coordinate condition becomes 
a fundamental law of nature changing the character of Einstein's general relativity 
theory into a theory of the gravitational field, valid only in inertial coordinate systems. 
Other (like Papapetrou 1951) who based their research on Fock’s work, do not go so 
far, but they also regard the coordinate condition (1.6) as essential for the deduction 


of the equations of motion. - 

———————————————— a —-—— m — —AÀ— À——Ó 
2 Some of the ideas presented here are due to Mr. Plebanski and myself; they will appear in fuller 

form in a book “Relativity and Motion" which we are writing. 


" 
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In our first paper (Einstein, Infeld and Hoffmann 1938) we used a coordinate 
condition different from the harmonic one; yet we obtained the same equations of 
motion later obtained by Papapetrou and Petrova. In a few subsequent papers (Ет- 
stein and Infeld 1949, Infeld 1953, 1954, Teisseyre 1954) we carefully analysed the 
problem. 

I maintain that the equations of motion have nothing whatever to do with harmonic 
coordinate conditions; they have much to do with the method of approximation. 
Indeed it is this method which uniquely determines the equations of motion up to 
the post-Newtonian order. I will show explicitly later that the violation of the harmonic 
coordinate condition does not change the post-Newtonian equations of motion as 
long as we stick to the approximation procedure. 


§ 2. Some notations and mathematical preliminaries 


We have a world line £*(t) and a field, say a scalar field y that depends on coordi- 
nates х^ and on time x? = t and also on the &*(#) and their time derivatives: 


p = (ай, t, é, £); EÈ = ——. (2.1) 


We assume that near the line £*(t) the field becomes singular and has the following 


form: 
е; 


ф aem. - puse 
gud "n T 9-9. (0 — &) + 9 Psr (а — e) — 6) +... (2.2) 
Неге 
e? = (x§— ë) (xt E) 2.3) 
and, therefore 
zi 
3 9 
ZUM м 
АРЕ 
Similarly ' | 
a1 
Ps = : go ts 
ET PT "S (2.5) 
We must distinguish between 
DE. - Әф fT БЕ. 
Ps and 9, = зв = Pe | (2.6) 


which generally are not equal to each other. 
If ф were not singular in the curve £*(t) we could have defined ф in the following 
way: 


Я =f Фд (X = E) daa. | <А (2.7) 
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Now we can narrow the definition of Dirac’s 6-functions in such a way that (2.7) 
remains true even if o has a singularity up to the Аһ order. Indeed we can show 
(see Appendix A) that such 6-functions can be constructed as limits of ordinary 
functions. This is convenient because, by the use of such ó-functions we can get rid 
of infinities without recourse to the renormalisation procedure. Thus all the 6-functions 
used here will have the property of changing into 9, where фіз a continuous function 
of the £*, £*, £ 

We can, therefore, use (2.7) as the definition of д, where д is the three-dimensional 
Dirac's 6-function, satisfying the following conditions: 


1) d(x) can be treated formally as a spherically symmetric function for which 
all the derivatives exist. 


2) d(x) = 0 for x 5 0. 


3) For every continuous f(x) in the arbitrary region Q(%,) forming a neighbour- 
hood of xg, we have: 


i dix O(% — xo) f(x) = f(x). 


Q(x) 


4) For an arbitrary neighbourhood (o) of the point xy = 0, we have 


Г dax 969 B[* = 0, for p =1, 2, ..., X. 
2(0) 
It is the fourth condition that distinguishes this 6-function from the usual Dirac’s 


6-function. The proof of the existence of such functions is given in Appendix A. 
(See also Infeld and Plebariski 1956 and 1957). 


Thus the *tweedling" means two things: firstly, the singularities are ignored; 
secondly, for x? the é* s are introduced. 


We can now return to the problem of the difference between g,, = s and 
p,- We have from (2.6) 


Va = Pat + 9. ; (2.8) 
This means: g,, and 9, are equal if and only if 9,5 — 0. This is certainly so if the part 
of ф that gives a contribution to y does not depend on &°. 


One more formula that will play an important role later and which follows from 
the definition (2.7) is: 


_ (2.9) 
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Let us now assume two functions e.g. p and y: 
g 


p= 24 + ASI В) P) s 


a 
- 


—3 —1 


yc ETE Gr P) (2.10) . 


We have in this case: 
Фр = фр (2.11) 


Indeed ee gy and ignoring singular expressions and those that become 
zero for 2° = &, we are left with the odd expressions of the order zero in (X? — £5): 


-1 as — — — © s)(xr — Ёт EN 


PYF Ф Y 
The parte of these expressions gives us: 


mA 1.87 eel а oe 
фф = f дез) diax E V + (pvd wv) sk 


0 


+5 imp Cae юе]. 
e? 


The last two expressions under the integral sign give zero, because they are product 
of symmetric 6-functions and odd powers of (x — &'). Thus we obtain: 


py = [$$ gy (зух = v. 
But (2.11) would not be true if ф or y had a singularity of order g~?. Thus (2.11) is 
true if о appears to an odd power in the singular parts of ф and v. Therefore we have 
to be cautious in applying the last equation. 
Thus we shall have vectors and tensors defined only а along the curve, like 


5°, T% etc. We can define the metric tensor along the curve Sap and g% and assuming 
(2.11) we have 


Sap 8 = Bu 8e = OF (2.12) 
To such tensors we can apply tensor algebra and tensor analysis but only along the 
curve. 
Since 


Qs) 
(where Q isa small three-dimensional neighbourhood surrounding the singularity), 
Sap will be a tensor along the curve if dis) (зух is an invariant. Thus we must discuss 
the transformation properties of дез). 
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The four-dimmensional relativistic Dirac’s ó-function доду is a scalar density 
because of the invariant equation 


] 5i dus = 1. (2.14) 
Now let us take 
дау = да) (x* — E°) (2.15) 
where 
&* = £^ (A) (2.16) 
À being an invariant parameter. Then we can form an invariant density function 
+оо +оо 
Јака = |р Es disse es дез) ` (2.17) 


As far the transformation properties are concerned this is the definition of Os. Thus 
дозу is the zero component of a density vector. From this definition it follows that: 


+оо 
ад dt 
fos d(g) x = IL» up Яя y= NET dax dt = fos dux. (2.18) 


(s) (3) AC) — со Qla) 


Thus the space integral of ô is an invariant. 
One more remark concerning the notation: If we have many curves we shall 


a . 
distinguish between them by the index written above: & ; a = 1, 2, ..., р. Then ме 
a 


should also write Ф, meaning: 


"EIE д (зух; д = дүз) (x — £5 (2.19) 
de | 
However, for the sake of simplicity we shall not write a above the “‘tweedle” always 


1 
understanding that it means the first (or only) curve: 9 = g. 
§ 3. The gravitational equations 
The gravitational equations expressed in contravariant tensor densities are: 


G^? = ра — = go? В = — 8л Т. (3.1) 
1 2 
In our case of, say, two particles moving along lines &* and & we have: 


ibi ama 
TË =t 0 Ht” ô. (3.2) 
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Before going further we shall introduce an invariant and finite line element, concerning 
the first particle (without writing the “опе” above). 


ds? = gy di^ dé. (3.3) 


From the tensor density we can now form a tensor along a curve (the first one) 


dt ~, dt 
р T ^? dix = 198 — 5 (3.4) 
ds i E ds 


Qla) 
What is "6 ? We shall see in the next section that the answer to this question 15 
a unique one; that the consistency of the field equations demands the following 
equations to be true (Tulczyjew 1957): 


108 = = u £v EP. Е = = £^; u = rest mass (3.5) 
s 


We shall show (§ 4) that not only has #8 the form (3.5) but also that the rest mass 
р is constant. For the moment, however, we shall use (3.5) as an assumption, without 
stipulating that м in (3.5) is a constant. 


From (3.5) follows: 


бий те; ВЕТ (3.6) 


and 


(3.7) 


The right hand side of the gravitational equations is uniquely determined by the 


condition that T°? depends linearly on the д". 


§ 4. The general equations of motion 


As a consequence of Bianchi’s identities we always have 


1 

Gp = БЕ) = — 8л т = 0 (4.1) 
3B 

where the semicolon indicates covariant differentiation. Taking the integral of 7, 

over the three-dimensional region surrounding the first singularity and multiplying 
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dt 
by =, we have: 
ds 


dt 2 
ds [ T% dox = A* = 0, (4.2) 
= | т", 
Nla) 


where 4% is a vector defined through (4.2) along the first curve. Generally we shall 
have 


a 


Ao 0053 со (4.3) 


if p is the number of singularities. Thus іп (4.3) we have as а consequence of the field 
equations, 4 p equations which we shall call the equations of motion of the p singular- 
ities. Indeed they contain 4p unknowns: 


P() and т; .s—1,2,8; a=1,2,...,p (4.4) 


Let us now write out 4“ explicitly (omitting the “опе” over the “tweedle”). We 
shall start with: 


т = T, + ra T (4.5) 


Introducing here (3.8) the first right hand expression becomes: 
jx eid 151 


тч = Т“, + T” o = (т), + (тебе > (46) 
+ similar expressions concerning other particles. 
We have: 
dt 1 1 1 1 dt 
c | (т ё #0), da) x = т =f. T^. dax = 0. (4.7) 


dta 
Al) 


This we can see even without explicite calculations by changing the volume integral 
1 


into a surface integral which must vanish because 6 vanishes on the surface of QQ. 
Therefore what remains of the integral of (4.6) is, because of (3.7), 


1 1 
rA a odayx = | т tox | 7 a, т т ёа) = (иё)? (4.8) 


Ms) 
20) 
Thus because of (4.5), (3.6) and the last equation our equations of motion (4.2) 
become: 


cM. 


us £o qe и Eo” + u [= gw ex =Q. (4.9) 


f 
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d Е 
We shall show now that w = SN 0, that is that the rest mass is constant. 


Indeed, let us assume, as we shall always do, that: 
а) g 10 
po —8 Гоо, o]. (4.10) 

Then, multiplying (4.9) by g,, £^, ме have: 

~w » , » эр — 1 —/ > , э 
Hg. ОЕ DE 5 № Ew, о ё gv £e — 0. (4.11) 
Because of (2.9) this is equal to 

Ов : 
M = Be BE) 05 р 0 (4.12) 


which changes one original equation (4.9) into 


(4.13) 


Thus we see that the rest mass must be constant. 

If not for the “tweedle” above the Christoffel symbol our equation (4.13) would 
be that of a geodetic line. For each particle we have four such equations of motion. 

Let us put back ¢ instead of s into (4.13) using (3.6). We then have: 


dm hos re. 
и ео 


~ 


AQ ё) т " auem moe (4.14) 


As we said before, we have 4 p equations determining the &^ and the m ‘s. From the 
last equation in (4.14) and (3.8) it follows that: 


m celsior о. 
F exp Да! 2а EP dt X (4.15) 
since 
m "os |a a: m 


From (4.15) we see that: 


t 


ds — dt IERI £g EF dt. 


о 
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This connection between ds and dt follows from the equations of motion. The normali- 
sation is such that for t = 0, both dt and ds are equal. 


Substituting the values (4.16) for in the last three equations (4.14) we obtain 


a 
the 3 p equations of motion for #°: 


a t ? PESO T Hoe | ый 
Feed pe 0; 2 = 1 (4.17) 


This equation, however, is not suitable for establishing a connection between 
the equations of motion and a variational principle. To do that let us rewrite the zero-th 


equation of (4.14) using (3.6) 


ато -— [24 io Ев 4.18 
(le) (eF), «pj* 5 4 
2 v 1 
от, A 22 о. 
ds 
ee РН | 4 
log Bas iin 3| = — GP Py (ean PEL 
5 20 
арек 
2 |. 1 2 (4.19) 


Therefore, the three equations (4.17) can also be written: 


«B 


Let us write, for short: 


£u k | а - (о £531 (87, Eo £6), V 0. (4.20) 


T- (eH? IP (4.21) 


Then we can rewrite (4.17) in the form 


£ — (log Ё), € + 12, Se Lc € (4.22) 
Indeed the “zero” equation (о = 0) gives the known equation (4.19). We shall multiply 
(4.22) by g,, and assume, from now on, that for 8а and their derivatives (2.17) is 
always valid, that is the ,,tweedled” product is equal to the products of the ,,tweedled 


expressions. Then we have because of (2.9) the three equations of motion: 


mE : = 1 dis 
(£c 20 = once (log С). ao 9 Sabor ВЕР == 0: (4.23) 
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This form suggests the existence of a Lagrangian: J£. We wish to see whether 
(4.23) is indeed equivalent to 


d ƏL ƏL (4.24) 
dt д IF 
We find: 
94^ ЯЕ 
з= Pi Rabb 
2.6 1 [ag 
ry (ее Bil TON 
E Е ge dne a E) L + See ICE) (4.25) 
dt 98% od 


Therefore the Lagrange equations (4.24) become: 
Ee TIR: pe l ae 
(Ear 6), о — Sar È (log L ),о— 5 Вараб 5^ + 
1 /8 ap я RA 
+1 (28:6 go Fs) Fo (4.26) 
2 \9é 20 
Comparing (4.26) and (4.23) we see that Lisa Lagrangian if 
Land ~~ 9g 
бов = BaB,k = T D 


28р 
d — 0. 4.27 
2g (4.27) 


: Els sf 1 
This means: = is a Lagrangian if the relevant part of g,, does not depend on or &. 


We shall see later that the first condition, that is the independence of the relevant 


1 
part of g,, from & is not always fulfilled. 
Thus (4.23) follows from a variational principle: 


д i Ёа — à E 4 [ y diy ду (2° — E) gap € 2] = 0, (4.28) 


t=—0o 

_ if we treat вов as function of x^ only that is if, while varying the function under the 
1 г 

integration sign, we ignore the possible dependence of Zap on & and &. Then, indeed, 

as can also be shown by a straightforward calculation (4.28) is equivalent to (4.23). 

If we call the result of such variation the ooo line" then (4.23) is the equation of 

a “geodetic line". 
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From what has been said it follows, that for a test particle — that is Фи > 0 
— 
and gag is not singular and does not depend on £ or £ — the equation (4.23) is that 
1 


? 


of a geodetic line and the “~” only means the substitution of & for x’. 


One can now show (Tulezyjew 1957) that the equations of motion follow from 
the field equations once we assume the linear dependence of 7° from the 6’s, that 
is one can show that the form (3.6) for t® and (3.7) for T isa consequence of the 
field equations. 

Indeed from Bianchi's identities follows: 


dt 

— Өт в d(ayx == 0, (4.29) 
ds 

1 

90) 


where 0 is an arbitrary function continuous on the worldline (2). Or, again omitting 
the “опе” over the ¢%’s and Ó's we can write: 


dt 
E 0(19,0(5));5 dax = 0. (4.30) 
$ 


1 
Aa) 
Thus this equation leads to the following one: 
A? 0 + А% 6,, = 0. (4.31) 
Let us start by calculating AŤ. Because 
à, —— 6, Ё . (4.32) 
we have: 
~ x dt DOT 
9,5 AP = == (— 19 + 190° )0,,. (4.33) 
ds я 
Since ¿0 = 1, this can be written: 


0,5 A76 = X cere ËB) б,в. T ^ (4.34) 
S А 


Because 0,5 is arbitrary, we have: 


EN 
& 
| 


dt Ч 
—— (— t98 + too EF), (4.35) 
ЕЯ ( F ) 


Putting here æ = 0 we find: 
28 — 190 EP, (4.36) 
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128 — 100 Ба ЕВ — m Еа ЕВ | (4.37) 


which is the proof of our theorem. Obviously A* = 0 gives the equations of motion. 


therefore generally 


§ 5. The approximation method 


We shall solve the field equations and formulate the equations of motion explicitly 
by means of the approximation method which we shall presently describe. 


Let us assume a function developed into a power series in the parameter 4 = — 
(for the moment we may assume с arbitrary and not equal to “опе”). 


0 1 2 


The indices written below indicate the order of А absorbed by the q’s. 

If the function g(x“) varies rapidly in space but slowly with x? then we are justified 
in not treating all its derivatives in the same manner. The derivatives with respect 
to x° will be of a higher order than the space derivatives. We can formalise this proce- 
dure by assuming: 


д 


— (9) =Po =Y, 5.2 
ma Pea a oe 


that is: differentiation with respect to x° raises the power by A, absorbed by the g s, 
by one. 

The problem is now with what order we should start the power development of 
the quantities appearing in the field equations? 

The quantity & will start (by an obvious convention) with the order “лего”. 
It will be an unknown quantity determined by the equations of motion; we shall for 


the moment not develop into a power series. Thus £* will be of the order- *one" and 


E of the order *two". 
From this follows: 


wn at ET es (5.3) 
4 6 


To begin with K is not pure convention. Indeed in the Nee approximation, 


which we hope to zd we have in the chosen units 


mass X mass 


mass X acceleration = —-—— — —- 
(distance) 


Since the acceleration is of order two, the orders of both sides will be equal, 
only if the order of mass is also *two". ^ 
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a 
In all the power developments we shall take into account only even (as in m) 


or only odd powers of A. (Infeld 1938). 


Thus because of the order with which we start m and A we have: 


700 — 700 НЕ то 3E 700 + BA 
2 6 


4 
ED c G p E G E G Е (5.4) 
3 5 7 


ee Ете 
6 


Now as to the g's we write: 
би» — "luy "s hs; g7 TE TM SOS (5.5) 
From the gravitational equation follows: 
l 


where 


Rag Rag у—== 


а RET В (t les " el "e By 


From the right hand side of equation (5.6) follows: Roo and R,,, start with the 
order “2” and R,, with the order “3”. The lowest order expressions on the left hand 


side are: 
; il 
for Roo : — pi Dos 
1 1 1 1 
for Rom: = pe E оня = E limes E 9 Е se 9 hio 
l 1 Л l E 
MIC c E > — > Путя: 8 
for T 2 (боен x 2 me 27 2 еа xi р) Боот 2 h x (5 ) 
Therefore 


hoo = hoo F hoo Foe 
2 4 

hom = hom + hom + ++ 
3 5 


be LT Л ит a Рп ар a (5.9) 
2 4 
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All the functions that will appear later are obtained from the A's by summation, 
multiplication, differentiation. Thus to every component, the following rule applies 
throughout: any component having an odd (even) number of zero suffixes will have 


only odd (even) powers of A in its expansion. 


§ 6. The Newtonian equations of motion 
We shall try to find the equations of motion in the lowest (Newtonian) approxi- 
mation. We shall do it in such a way as to make the generalisation to the post-Newtonian 


approximation as simple as possible. 
Because of (5.6) and (5.8), the field equations of the lowest order are in Ag: 


P 


1 аа 
== == hoo,ss = — O07 (тт) то mo, (6.1) 
2 2 2 2 2 2 а=1 
ог: 
р аа 
hoo, = 8% У) тб. (6.2) 
2 a-12 
As the solution of this equation we take the Newtonian field, that is 
11 22 
hgg— —2mr'—2mr3—... (6.3) 
2 2 2 
where 
r? = (x! — &) (x! — Ё). (6.4) 
We shall also write briefly: 
hoo = Ф (6.5) 
2 
and in the two body case, which we shall assume for the sake of simplicity? 
Г 
т 22 
= —2тг 1; ——2mr! (6.6) 
1 2 
= БЕТ 
It is also worth remarking that because of (5.5) and 
g” Epo = д", 27 (6209) 
we have 
h — —ф (6.8) 


2 
Е EE St a УЧР ЧИДЕР А or A c Г 
3 There is, I hope, no danger in confusing this g in (6.6) with & = [Вов 


On 
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We see that the relevant part of g, that is g does not depend ор £; therefore we 
have: 


— — ~ 
оо» = hoops $»; 
2 2 


32 | 
9 86 (6.9) 


There was по reason to assume that m is а constant. This follows, however, from (4. 15) 


or the first equation in (4.14). 

Because of (6.9) and because in the lowest order we do not have any products 
of the A's, the conditions (4.10) and (4.27) are satisfied, this means that the path is 
а "geodetic line" and can be deduced from a Lagrangian: 


hc ПГ Ne i ETE. 
L = (Bae &* £^); авв) (6.10) 
E. 
1 2 = т 
О ИИ (6.11) 
| т) 
Thus the equations of motion, up to the second order, are: 
dob PL eg. (6.12) 
dt IES DES 
that is: 
r 2 
x m 
yee —]. 6.13 
ENS ue 
Past . . 
The Lagrangian £** for both particles is 
DES ; 
Im. 29525 
Т = 1 тё тарар T (6.14)) 


The Newtonian Lagrangian for both particles is of the fourth order! One more remark 


dt 
Since a = A , we have here: 
m=u 
ИРА а (6.15) 
ме тацда | 
Е r 
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We can also find the equations of motion directly by going back to (4.17). We then 
have simply: 


+l i} qo (6.16) 
оо 
Or 

t 2 E 

ОЕ е (6.17) 


which is identical with (6.13). 


§ 7. Transition to the next approximation 


To find the equations of motion up to the fourth order besides hoo we must also 


2 
known 


h 


mn ? 


Ie as (7.1) 


от ? 
Thus we must proceed to calculate these expressions. As far as the first two are concern- 
ed, the task is easy. The left hand side of the corresponding equations is written out 
in (5.8) and the right hand side is given by (5.6) and it is: 


13112 2 
for mn: —4лд„ (тд + m д) 
2 2 
ПЕ he be ear) 
for om: 8л (m £" ô + m n” ô) (7.2) 
2 2 


Therefore, for À,, we have the equations: 
2 


1 1 1 1 1 l 
ESE 9 Рт zs 9 Ата xis 2 Pens ms T DI him =f о; smn = — ә б.н Ф, + (7.3) 


We аге looking for a solution of a Newtonian character. Such a solution is: 
Pern = дни Ф. (7.4) 


The choice of these Newtonian solutions for hoo and Ay must be regarded as part 
2 2 
of our approximation procedure. More about this later in § 9. 


The next step is to calculate №„. Again (5.8) and (7.2) gives: 
4 


о (aves PURO 
2 e р) oss F а = ол (тё ô + ия 6) (7.5) 
Introducing (7.4) into the last equation we have: 


1 h 1 h Nee | 22592 . 
Hire 2 оњ Hy Позз — Pono = 8л (тё + тад). (7.6) 
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The solution which we shall adopt here is 
—2ft"—2gm'. (7.7) 


This is not the only possible solution. About a more general solution and its influence 
upon the equations of motion we shall say something in 89. 
The calculation of hog is more troublesome. We shall give it in Appendix B. 
7 


Here we shall quote the relevant expressions of Лоо, that is those that give а contri- 
4 


. = ned ry 
bution to hoo and hova They -are: 


hap ~ 24m) ar — 37 p "RO Tee Seow (7.8) 


Now we should like to see whether these expressions for A 


p Лол. liop are пер 
3 


4 
that they make the equations of motion equivalent to those of a ,,geodetic line", that 


is, if the conditions 


TARTE (7.9) 


cron О (7.10) 


are satisfied. Both conditions are satisfied for h,, and й„. Indeed this expressions 
: 


have a singularity of the order L therefore an odd singularity, therefore for them 
r 


(7.9) is satisfied. Furthermore in their relevant part neither € nor & appears expli- 


citly; therefore for them (7.10) is satisfied. But this is not true of Ро. It has an even 
4 


1-2 
singularity because of the appearance of @) . Yet this does not matter since Aq, 
K 


appears only linearly in the equations of motion. Therefore we can disregard the 
condition (7.9) for ho. But is the condition (7.10) satisfied for hog? It is certainly sati- 
А 


DEN (2) мо (5) 
sfied for the first two expressions, that is for 2 Col rj — 3 р т т), ѕіпсе nei- 
ther of these two expressions depends on &, &. We must be more careful with the 


third expression: 
2'2 212 a 21249927 
— MT,99 = — MT, ps nn —mram. (7.11) 


The first еа on the right hand side does not черо on Ё, È, but the second 
1 
expression contains ү. Here we may introduce for т р = — т E, the Newtonian 


Li 
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value since the mistake will be of order "6". Thus we can put: 
QUO à AD 12 2x -1 
Qe e T TEE тга Е = [* m (x3 — 75) (7) ] (a5 — n)r-*. (0.12) 
Г 9&5 


This expression, differentiated with respect to x” gives zero at the point x! = &. 


Thus 
a, = 0 (7.13) 
But 


A similar situation occurs in the last expression in hog: 
4 


AyD 
<2 nem G ЭН (7.14) 
We have: 
A 1 2 dps 1 2 
B,—2mmr[(r], =2m mr Цит. (7.15) 
- 12 
4 = 2mm (г), 
Therefore: 


В, = 2 B. (7.16) 


We can, however, easily find an auxiliary field which we shall denote by hoo 
such, that: 


2. = 
oom a Боо, (7.17) 
From this it follows that such he, is: 
4 
ae 2 1 2 
hog = 2 (т)? г г ЗИ mr т, gen n Е тт (7.18) 


4 


$8. The post Newtonian equations of motion 
If we look at the general equations of motion (4.17) we see that — up to the fourth 
order, hoo will appear only once, that is in the expression S This means, because 
4 


of (7.17) that a Lagrangian up to the fourth order exists*: 


L* = —— (8.1) 


4 Compare Fichtenholz (1954) where the Lagrangian is found mechanically from the explicite 
equations of motion and not as here, the equation of motion from the Lagrangian. 
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with 
ТА 
| di ) ge is 
in which only g% % Zoo and 
goo = Noo + hoo + hog = НФ 4- hoo (8,3) 
2 4 4 
Therefore: 
T А, ш, M t ni бу 
L* = (Lp + hoé обе, Е). (8.4) 
3 
Here we introduce the values: 
Ss 2 m 2 E 
qe 2 mrst: А =4mr Ni 
3 
= 2 dE 2 IEA 1392 
hoo = 2 (т)2 172 — 3 n n mri — m rse TA +H mmr. (8.5) 


Remembering that up to the fourth order 


i 1 1 
1 +а-+ а)" =1 - : 
(+a+a*=14+ 5+a— ga (8.6) 
we obtain from (8.4) and (8.5) the Lagrangian: 


2 


~y ~ ~ т 2 . 5e 
L7 = ES utem ББ тті 5 тг (Es & + ns 28) + 


Nj = 


2 1 2 
e MEAS Borm 1 т (т - т) Lie? HA 
-1 ЗЕНА $ E52 aad ee cadi r nS nr 
о z ^nse TU. (8.7) 
This is supposed to be the Lagrangian for the first particle. But we wish to find 
the Lagrangian for both particles. This means: first: a Lagrangian which gives the 
same equations of motion as the Lagrangian (8.7); secondly: a Lagrangian which is 
invariant with respect to a transformation, changing: 


1 2 : 
mé 2. mm. ' (8.8) 


327; 1 
Let us multiply .£* in (8.7) by m. Then the equations of motion will, of course, 
be the same. Let us add: 


1 orate eae A ERN 
—3,mmwv-—gmGtmy. (8.9) 


Then, with this addition, obviously, the equations of motion for the first particle will 
still be the same since there is no contribution from (8.9). Now with these changes, 


d 
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the only expression not invariant with the respect to the change (8.8) will be the last 
one in (8.7), that is: 


1 2 
== Dp m mM Го т 7’ (8.10) 
But instead of it we can write: 
| 12 : 
> mM M TyEsnr Esn / (8.11) 


which is invariant with respect to the change (8.8), and gives the same contributions 
to the equations of motion. 
Indeed from (8.10) we have the following addition to the equations of motion 
of the first particle: 
Ter 1 349 
z mmr, gem р n rs m mr, gare р 5. (8.12) 
From (8.11) we have the addition: 


Th 172 1 1—2 
y тт (ти "c 9 m mr, gap gn Bay P 


(ouem T ss 
F m m (r, gares NESET, ennas NNT (8.13) 


S ee. Tug муб 7 
- г, гп ур HI" — Г, es gr en ES 7’) = 5 m m (Г, gn ит Ее) 


But the last expression in (8.13) equals: 


rog 3 te tee $ (En — ар) (& — ар) 
И Т Тут yr 7 = Sues ne L r3 r (8.14) 
Since 
Sd m GU r 
ТЕА аи 3 т) (8.15) 


we see that (8.14) vanishes and, therefore, (8.10) and (8.11) give the same contribu- 
tions to the equations of motion. 


Therefore, if we call the final Lagrangian for two particles Le, we have: 


аль Сън опат mm 
2 2 r 
1 2 н, 2 E 
3 m m 


-3 (Bé п) +4" 


s Telus na. 
T isy ее (8.16) 


12 
1 ти (в m) TE 


I.: 
a eee т)? + TUTTI + z m M T, ва Es ур 
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By a straightforwards calculation we find from this Lagrangian the equations 
of motion of the first particle: 


2 1 
z e 2 d ONE Pee: 
LIE =m ЕСЕ д © | 4- 
г | ап 2 r r P lagz 


-+ [AES (0 — ® 341] (2) Qo улет. (0.17) 


The equations of motion for the other particle are obtained by replacing 


1 2 2 8931 
Mit, Ce by mcm nee (8.18) 


respectively. 

The generalisation of this result to p particles is almost trivial, if we take into 
account the changes caused by the addition of these particles in hee and outlined at 
the end of the Appendix B. These additional expressions are due, say, in the case of 
three particles, to the interaction between the second and third particles; that is in 


the equations of motion for the first particle they will give a contribution proportional 
LEE 
to т m m. These expressions appear in .É** from two sources: from А7 and from ф g? 
4 


in (8.6). If we now denote the “distance“ from the a-th to the b-th particle by: 


(ab) a 


b a b Б 
r= [(£— 2) (E — 8^ ИК 


then we have the Lagrangian for the р particles: 


д РИН 1 bs LY Se pe 
мары т vm 
aa a a b 
Y NE Е E 4 Е à) nr XS Y ae 9), fis 
dni аа 1 
(ab) 
AIC pira NIE yes aa ste 
тл 
E. у са b iy ee 5 5 " x 5 i g 2 Са пси pec Ак 
e bir (8.20) 
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In the case of two particles the Lagrangian (8.20) reduces to (8.16). The only new 
expression appearing in (8.20) is the last one; in the case of three particles it is equal to: 


КОЛЕЯ п 


if we accept the notation of Appendix B where бод is the change in hy, caused by 
4 4 
the interaction of the second and third particle and / is for the third particle what f and 


gare for the first and second. Thus — : р kis the contribution to the Lagrangian of 
the interaction between the second and third particles coming from — s g? in (8.6). 


8 9. On the choice of the coordinate-system 


‘The harmonic coordinate condition is 


Vetere e acsi (9.1) 


which in our case means 


M" ess аа Оле 0 (9.2) 
2 3 2 
None of this conditions are fulfilled in our coordinate-system. 
The values for hog and A,,, accepted by us here were the Newtonian values. Our 
2 2 


convention is that their choice characterised our approximation procedure. Yet, with 
some justification, this approach may be regarded as too formal. Indeed, instead of 


our values for À,, we could have chosen 
2 


» = af 
h ИЕ: hinn Е Ч тп ae Gaim 
2 2 2 2 


the a’s being arbitrary functions. This change from h,,,, to X „n could also be induced 
2 2 
by a change in a coordinate-system that does not disturb our approximation procedure. 


What is the physical meaning of the choice а, = 0? It means that we assume the - 
existence of a coordinate-system in which each of the two bodies reveals its spherical 
symmetry; a coordinate-system in which for m - 0, Е — 0, the field goes over into 
that defined by the Schwarzschield solution in an isotropic coordinate-system. The 
choice of such a coordinate-system is implicitly assumed by our approximation proce- 


dure. However this choice of a,, = 0 refers only to the beginning of our approximation 
2 


procedure; therefore it does not refer to Ао. Indeed if we replace hom by 
3 3 
Pom E hom Е Lorm (9.3) 
3 3 


ао being an arbitrary function of x*, then the equation (7.5) is fulfilled just as well. 
рр 
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This change in Ag, сап be induced by a change in a coordinate-system which does 
3 


not disturb our approximation procedure. Such a change also induces a simple change 
in Aoo (see Appendix В) 
4 
) | $ 
П? оо = hoo + 2 ао (9.4) 
4 4 


Therefore it would seem that the Lagrangian and with it the equations of motion would 
change. The expressions that change in the Lagrangian (8.4) are: 


— 54 
hasë + > hoo (9.5) 


and the change induced by them in the Lagrangian is, because of (9.3) and (9.4): 


о ae ae 
ИР а ae (9.6) 
Therefore 
5 [4 £* dt = ба) |= 0. (9.7) 
ti ti 


This means: the equations of motion are uniquely determined up to the fourth 
approximation by the field equations and by our approximation procedure. Neither 
the harmonic coordinate condition nor any other coordinate condition played any 
role in our derivation of the equation of motion. 


§ 10. The general theory 


Let us now formulate the general theory®, according to which we proceeded in 
our special case, in which we found the equations of motion of the second and fourth 
order. Of course, such a general theory is of little practical value since there would 
be hardly any physical meaning in developing the calculations one step further. 
Moreover, it seems — and we shall discuss this later —that by the proper choice of 
the coordinate-system we can annihilate all contributions to the equations of motion 
beyond the fourth order. Yet from the formal point of view it is important to know 
that the procedure can be pushed as far'as we wish. Of course we do not know anything 
about its convergence. 3 

Before we formulate the general theory let us recall what has been done is 
We had the Newtonian equations of motion: 


СА IS » dig x = A US E At: d(gx = 4" - = 0. (10.1) 
ds ds 


5 The ideas presented here are a few years old. Scheidegger (1953) refers to them in 8 5 stating that 
they were suggested by me. I found a more explicit formulation of similar ideas in a thesis by Rameswararao 


(1955). 
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1 1 
Since m appears as a factor we called these equations (after dividing them by m) 
1 2 


the equations of the second order. But here, since it is multiplied by m it appears as 
2 


an equation of the fourth order. For this section therefore, let us rename the order 
of the equations of motion calling the Newtonian equations of motion those of the 
fourth order and the post-Newtonian equations those of the sixth order. Let us also 
put generally: 


dt ~ di 
= T. d ЖЕ Ao Mer 25:0] x= £. P 
ds | 291^ 5) 2—1 ds М x (10.2) 


Thus for our post-Newtonian equations of motion we have: 
AS 354 i Qi CAECIS () (10.3) 
3 5 2 4 
These equations gave us: 

a a a a a 
тет m сс en (10.4) 

9 MP 0 

a 


where &° is the motion in the Newtonian approximation. To find these equations expli- 
0 A a a 

citly we used the Newtonian equation of motion in A” since the use of & instead of & 
6 0 

would give a contribution of the 8-th order to the equations of motion. We express 


this idea in symbols and write instead of (10.3): 
A isn 5) РЯ A (©) =). 
4" Gt + =") аР 4" G = ©. (10.5) 


Thus Ange -+ 8) also gives a contribution of the order „six“. The field was solved 
so as s dues (10.4) in: 


m, n up to the order A,, (£) 
2 0 
о, n 5 33 >> 25 hon (é) 
3 0 
70. о 22 2 2 2 bs (2) 
4 0 


Let us assume now that we wish to go опе step further. We then have the equations 
of motion: 


3 


AV Et bye Ao fet E OU 0, 
0 2 4 5 0 2 T 0 


Á" (E -- £ 4- £) "Е" © — 0. (10.6) 
4 0 2 4 6 0 2 8 0 
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In AS, AS, А", А", the argument in (10.6) is different from that іп (10.5) therefore 
Sheree БКС 


they give contributions up to the eighth order. But to find A"(£) we have to know 
0 


hon hoo» (10.7) 


6 


h 


mn? 


all functions of é. Thus simply denoting: 
0 


9% = — 8л (rs — 5 g^ r) (10.8) 
we have to solve the equations (omitting the a’ s above the é, s): 
Re (E + E) + AE) = Q"" (E + £) + Q™ (5) 
OG Xe о Ac 2 4 o 
В" (E + в” (£) = Q™ (E+ & + Q (6) (10.9) 
Soa 5 0 q 0 X 5 0 
ЕО == Е ао (220 ою (5 


2 0 2 4 4 0 2 6 0 2 0 2 4 4 0 2 6 0 


This would seem to be an impossible task, since to solve the last equation (10.9) 
we would have to know &, which we wish to find by equation (10.6). However this 
4 
is not so, because we have 
il p aaa 
Н = — > hops = — 4n У тд = Q™ . (10.10) 
2 2 2 €: 2 


for arbitraty motion. Thus we may rewrite the last equation (10.9). 
R (E £)-4- КО (£) = QUE + £) + Q™ (5. (10.11) 
A 0 2 во До 2 e) 


Collecting here the expressions of the sixth order we find hoo(é). Similarly we can 
E- 6 o0 
find ho,, and h,,,. Thus we can push the approximation one step further. 
5 4 


We can now formulate the general theory. 
In the development of g,,; R^; T”, we took into account arbitrary motion. 
Under this assumption we developed, say 


Tm — Tmn a” T" E (pt ue (10.12) 
4 6 8 


But, instead of arbitrary motion, let us put into the arguments certain motion developed 
into power series: 


= ЕНЕРУ. gu 
0 2 4 
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and write, say, 
T(E +§ +8). re 


2s 0 2 
If developed properly, the above expression gives us contributions of the order 2s + 2 
and 25 + 4. Remembering this, let us assume that we have solved the equations of 


motion of the order 2r: 


ао 
DEED 


Ao (EE ER 

3 0 2 2r—4 2r — 6 2r-1 0 . 
о S 
4 0 2r—4 6 0 2r—6 2r 0 


This means, that we have solved the field equations in: 


m, n up to the order А„„ 


2r—4 

о, т y 2 22 2 h Om 
2r—3 

0, 0 39.195 2? >> hoo 
2r—2 


Now we wish to solve the equation of motion of order 2r + 2: 


4 Parva ie .+ A (B =0 


2r — 2r+1 0 
PLN CE LETT 1) (10.16) 
4 0 2r—2 2r+2 0 


The arguments in А, er E А”, nits А" ате different in (10.16) from those in 
3 2r—1 4 2r 
(10.15). Therefore they give contributions of order 2r + 1 and 2r + 2. But to find 


А” (&) we have to know: 
2+2 0 

Pann» ho, hoo: 
2r—2 2r—1 2r 


We have, therefore, to solve the equations: 


T ++. . + g"" (5 = поне гордое ЕО: 
2—2 0 


2r—2 0 2r— 


ger Е)" = em us E) +... +9": 
2r— 2r-1 0 2r—4 2r—1 0 
rE ај Budd. ome. 25 id put SF uc. 
r— 2r—2 
Everywhere in these equations, with the exception of R% and Q% we substitute 
2 2 


the motion already known. However R° = Q for an arbitrary motion. Then we 
2 


(10.17) 


can replace the last equations in (10.17) by 
nee та 4 ре ааш рт TH 24 2: um nid 


2r—4 
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In f, the expression Й„„ appears, for the first time. In R” the expression hy,, 
r— 2r—2 2r—1 2r —1 
appears and finally in R” the expression ho, appears. Collecting all the contributions 
2r 2r 
of the highest order in these equations and putting them equal to zero we can find 
Вп > Rom» Лоо; for hoo the equation is purely a Poisson equation! 
2r—2 2—1 2r 2r 
Looking back at equations (5.8), we see that if Aom h,,, is a solution of (10.17), then 


2r—1 2r 
» 
= | 
hom = hom T Loom 
2r—1 2r—1 2r —1 
, ди 2 | Д 

homn - hnn T п x ап" (10.18) 
2r 2r 2r 2r 


is a solution also. 
For example let us put r — 2; that is 


, 
h nr hom 4 Lom 
3 3 3 


3 
h mn h nà xi ann Е Anm 
4 4 4 4 


The choice of these functions can always be achieved by a coordinate-transfor- 
mation from a coordinate-system in which the a’s equal zero. Then in the equations 
of motion of the 8-th order, the derivatives of these four functions will appear. Generally 
they can be so chosen as to anihilate the expressions of the 8-th order in the equations 
of motion. But it is difficult to judge whether such a coordinate-system would have 
any physical meaning. In any case, up to the post Newtonian approximation, the 
choice of the coordinate-system does not play any role as long as we stick to our approxi- 
mation procedure by which its beginning is determined, that is the choice of hoo 

2 


and h,,. 
2 


APPENDIX Ав 
To distinguish between our 6-function and Dirac's 6-function we shall here denote 
the former by ô. Our aim is to give a “realistic model, showing how to construct 


a sequence of (е) во that 6 = lim (е) and such that for every ё: д 
р £0 


fôle dar=l; | 8 (0) г? зк 0; p—L2..* (A.1) 


Such a model can be gained from a model (г) of an ordinary Dirac 6-function 


satisfying the following conditions: 


90-99-44 (2) | (A.2) 


€ 


в This is an abbreviated and changed version of two papers written in collaboration with J. Plebanski 
appearing in Bull. Acad. Pol. (1956 and 1957) 


# 
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where 4 = is such that 
E 
ET DP Јав z АА РОМЕ Leh (А.З) 
4л 
0 
always exists, and 
DO = | (e) dx = 4n | 2 A(z) dz = 1. (A.4) 
ò 


k 
If ó(e) does not have this property, it can be made to have it by multiplying it by (=) 
€ 


and renormalizing?. Thus with such 6’s we can form the model of our 6 function in 


k 
9 (e, r) = a (2) 70) (A.5) 


We have to show, that such a choice of 6, satisfies (A.1). To do so, let us start with 
the first equation (A.1): 


the following way: 


^ 4л [ 2 \* 
? k 2 —Dpor- 
зах = Т B 2 fae: dz = D it (A.6) 
0 


Now as to the second equation (A.1): 


k d k 
Ја i о ar [лова s (i ek- DG) = 0 (АЛ) 


0 
for every integer p, if 
l <p <k. 
Thus the s defined by (A.5) satisfy (A.1). 
This procedure can easily be generalized. We can introduce modified Dirac’s 


A 
functions ô from conditions 


Горе a (A.8) 


m а. ее N Го А. 
7 e.g.: Let us take д (=) = (270) lh: 2-3 ele! 57, that is A(z) = (27): e^: 2°, Such A(z) shall be 


k+3 -1 
changed into A(z) = (22)! 2 =) | ее zk ет» z* 
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where wœ, are arbitrarily prescribed numbers. Indeed the realistic г) satisfying (A.8) 
in the limit е + 0 is the following: 


A k 5 
дә = oa 3) - Не] (А.9) 
s=0 


We fina. 
обо feos Н 
fie +e = У ре. oy nano (SY tre frage 
=0 0 
k 
EU OCA PCR Vau 
= У (D&-5) ne 53 ek—? Г). (A.10) 

Sx 2 


This is different from zero and finite for € = 0 only for s = Е — p. We have: 
f 99 rds = ap + 0 (9 


and for € — 0 we have (A.8). 

We see that the use of Dirac's functions requires the prescription for the values 
of wp) s. The one used in this paper is the most convenient for our purpose; it requi- 
res Wy =0; p=1,2,...,k. 


APPENDIX B 
Коо up to the fourth approximation equals: 


1 1 


3 i 1 
len a oH =o ооп F Mos» os — 2 Фоо т og Pis Pos — y ФФ. (B.1) 


Therefore — 
RY = Rog (1 + 24%) = Roo (1 — 29) 


RM — RO V—5 = (1— 9) R= Rog (1 —39) | (B.2) 
In (B.2) exceptionally g = | Sap |. We have: 
i 1 
R” = — 9 Poss + ФФ, n peo : 9,00 — 3 ws. (B.3) 


The right hand side of our gravitational equations 


1 
RE = — 8л (ч — 3 g^? г) (В.4) 
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is up to the fourth order for the "zero-zero" component: 
| 1 1 1 | | 
m рее ТЕ о (4,00 grt. 00 nes В.5 
& [т z Oeo) Te Tos DM NE 


= — An (т% + 7% + T 
2 4 4 


Because 

12 
1 1:526 Um 
m—zmbb mm 
4 Ри тл Г 

12 
2 TRST "n 
ру 
4 2 з х 


we have for the right hand side: 


| 12 1? 
w a a oe 2595 1 5 
— 4n (mo imi пот (B.6) 
2 2 2 "SE 2 ГА r 
Now on both sides we shall take only expressions a) of the order four, В) those that 
give a finite contribution to А оо,„. Thus this relevant part of Aog satisfies the equation: 
4 l 


И: 
а +в (8 ттд д JI (B.7) 


M 


Because = | «<a > 


А вы ват; = о HOMES вв). 


we have for the contribution of C to hy coming from the ^ rst two ‘expressions: - 3 
A: > < RV < eri: у ra x "14 e ? LE юс - ~ ewm ER t mm e 194% x 
rác 2E 2 4 EN " 


(ur iow dol Esth 8B Ў аа. Bee Es. E 
ee iU elg а M MEME. cd 
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Thus in our case: 


12 
4mm 
Cis 3 (B.13) 
1/27 8 
Therefore: 
гаа 2 rM 
á ae: mm 2 
Поов = 800 T 6: £s + 8л s О 7 ? (B.14) 


Finally therefore, we have: 


hog >— т ое. (б): Gi — 3 dr р (QE MENS Hey (B.15) 


4 


Let us now generalize А оо for three particles, again only looking for expressions which 
4 


— 
give a contribution to Agom The only non trivial expression of this kind are those 
4 


23 
proportional to m m. 
(a b) 
We denote by г, the "distance" between the a and b particle?: 
(iy a b a b 
г) = (E — &) (&—&) (В.16) 
and we ask: what are the contributions to (В.4) coming from the third particle and 
20 
proportional to m m? We have now:: 
ф = Ло =Л-+Е-А . (В.17) 
2 


ая TANEDE m7 3) 
НИТ | i misc oer) ; facem EM ER 
Then the additional E for which we look in (B.3) are: 


1 f 
iG k),s Ec 2 S р et = kgn V 3 Soos (B.18) 


where So, denote the additional expression in Wn The additional expressions in (B.5) 
because of (B.6) are: 


— M mm ue Md (5 ate 5). (В.19) 


Therefore the additional expression to hoos is: 


и ею а 


From this we find 
S oS nda G dig addi € Y cen Circo вар 


a2) 
8 Previously r =r 
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1 
Therefore, as we see, Sy does not depend explicitly оп é. Therefore: 


S00,» = Soo, 7 (B.22) 
and 
RT 1 p я 3} (23)(31) E ay (12)(13) 
ттт (9) + (тг) (7 г) ] ea 


The last very simple question with which we shall deal here is the change 


from he, to т 
3 3 


h om == ho, sik Com: (B.24) 
3 3 3% 


Putting this into (B.1), we have 
h’ оо = hoo + 240,0 (B.25) 
4 4 4 


which is identical with (9.4). 


КРАТКОЕ СОДЕРЖАНИЕ 


Л. Инфельд, Уравнения движения в одщей теории относительноети и принцип 
наименьтего действия. 


В настоящей работе обсуждён способ формулирования уравнений движения 
для p тел, принимая: 1) уравнеия поля общей теории относительности, 2) факт, 
что тензор эиергии скорости отличается от нуля единственно вдоль р мировой 
линии известной формы, иначе говоря, что тензор этот линейно зависит от 
б-функции Дирака. 

Показано, что на основании этих простых предпосылок можно получить надле- 
жащий лагрангиан, избегая при этом обременительных исчислений. Лагрангиан 
этот ведет к вненьютоновским уравнениям движения. 

Выяснена тоже связь между полученными уравнениями движения и урав- 
нениями „геодезических линий”. 


REFERENCES 


Einstein А., Grommer J., Sitz. Berl. Akad. Wiss. 2, (1927). 
Einstein A., Infeld L., and Hoffmann B., Ann. Math, 39, 66 (1938). 
Einstein A., Infeld L., Ann. Math., Can. J. Math., 1, 209 (1949). 
Fichtenhole І. С. Journ. Exsp. Teor. Fis. 27, 563 (1954). 
Fock W. A., J. Phys. (Moscow) 1, 81 (1939) 
Infeld L., Phys. Rev., 53, 836 (1938), Can. J. Math., 5, 17 (1953), Acta Phys. Pol. 13, 205 
(1954). S 
Infeld L. and Plebanski J., Bull. Acad. Pol. (in print). 
Papapetrou A., Proc. Phys. Soc. (London), 64, 57 (1951). 
Petrova, Izh. Eksper. Teor. Fiz. 19, 989 (1949). 
Rameswararao B., Thesis,of the Banaras Hindu University (1955). 
Scheidegger, A. E. Rev. Mod. Phys. 25, 451 (1953). 
. Teisseyre R., Acta Physica Polon. 13, 47 (1954). 
Tulezyjew W., Bull. Acad. Pol. Cl. III 5, 219 (1951). 


Vol. XVI (1957) ACTA PHYSICA POLONIGA Fasc. 3 


ON THE FREE-ELECTRON THEORY OF ABSORPTION:SPECTRA 
OF SOME LINEAR CONJUGATED SYSTEMS 


By S. OLSZEWSKI 
Electrode Processes Laboratory, Institute of Physical Chemistry, Polish Academy, of Sciences, Warsaw 
(Received December 5, 1956) 


The AS FEMO method was applied in this work to calculate the absorption maxima 
in the spectrum of symmetrical cyanine dyes and polyenes. The results were discussed 
in comparison with experimental data and these of the simple free-electron model. 

In the case of cyanines considerable improvement of the theoretical values is obtained 
by introducing the resonance barrier in the FE model. The height of this barrier is deter- 
mined by the stabilization energies of the cyanine nuclei. The effect of the energy states 
obtained by the introduction of the barrier to the simple FE model is discussed in detail 
and the correspondence between the results of resonance theory for cyanines and those 
of a FE model with a barrier was indicated. 


$ 1. Some Results of the Application of the AS FEMO Method 


The recently developed AS FEMO (antysymmetrized free-electron molecular 
orbitals) method (Olszewski 1955) allows improved calculations to be made on the 
absorption spectra of non-branched bond systems of conjugated molecules in the free- 
-electron scheme. In this method even long chains may easily be investigated. 

Let us assume, as usual, that in the free-electron model, for the above molecules, 
the energies corresponding to the electronic eigenstates are: ; 


h?n? 
E = En = 8ml? X (1) 
where n is an integer, 
and the eigenfunctions are: 
y, = (2/L)% sin = S: (2) 


where L is the free-electron path. 
This may be called a first approximation solution to the problem. (The interaction 


potential of electrons with the core is taken as equal to zero). In the second approxi- 
i (211) 
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mation (which corresponds to the AS MO treatment) we built the eigenfunctions from 
the antisymmetrized products of molecular spin-orbitals, the mutual electronic inter- 
action operator being added to the Hamiltonian. We then can use the formula given 
below (Roothaan 1951) for the energy transition from the occupied ground state 
level i to the excited state level a: 


E(* Ф,) — E( Фо) = Н, + >, QJ, — Kj) — Hi 
j 


mits 2i (2J; Ki) | (Jia Ki) ir Ky, (3) 
J 


where the summation runs over all orbitals of the ground state; i.e. for an even number 
of z-electrons the summation is taken over N/2. The signs + and — correspond to 
the singlet and triplet states, respectively. 

For the elements appearing in (3) in the case of the non-branched free-electron 
path, the AS FEMO method gives the following expressions: 


Houck 
e? m? d n? 
Jam = = |= + (n3 — тп Si(Znz) s (n3 — m3)m Sina) A 
oie: 3 Si (2пл) | dese 
Jon = S P ad 9 n "y Kan, (4) 


Km S sitio ma (i mbm) + Sim — maa ( : em 


n--m nm п — т nm 


where n and m may obviously be exchanged. 


The scheme given above has been applied to the polyenes and symmetrical 
cyanines with various №. The wave lengths corresponding to the lowest electronic 
transitions (N/2 + №2 + 1) are calculated! in Tables I & II. 

For L the (N + 0:4 value was taken for polyenes, and for symmetrical суапіпеѕ — Nd, 
where d = 1.40 A is the mean length of the C — C bond in the chains. 

Theoretical results compared with the experimental data for polyenes show that 
_ the free-electron model cannot be applied to the treatment of these molecules. To Anax 
there correspond much shorter waves than those obtained by the theory, and the abso- 
lute value of the difference between Ags pemo and А.р: increases with increasing chain 
length. This situation is caused mainly by the well-known fact, that the z-electrons 
in polyenes do not form a complete free-electron gas, but are localized on the double 


SS as ЖЕ ЕВА a E RI. 


1 [n the case of | polyenes with N = 2 (ethylene) and № = 4 (butadiene) the exact data are given 
by Olszewski (1955). : 
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TABLE I 


The locations of absorption maxima for polyens 


1 (A) ГИ eA) 
Polyen N AS FEMO simple FE exp. 
| singul. triplet model 
Hexatrien 6 3170 4590 4530 2600 
Oktatetraen 8 3930 5800 5820 3020 
Axerophtene 10 4750 7150 7110 3460 
Anhydro-Vitamin A 12 5450 8220 `8400 3690 
Dimethylpolyen 14 6420 9930 9700 3995* 
Dimethylpolyen 16 7110 10990 10990 4090* 
Dimethylpolyen 18 7920 12310 12280 4160* 
Dimethylpolyen 20 8700 13570 13570 4230* 
В-Сагоѓепе 22 9500 14890 14870 4510 
Dehydro В-Сагоѓепе 24 10290 16160 16160 4750 
26 11090 17470 17450 
28 11870 18740 18740 
Dehydrolycopene 30 12670 20060 20040 5040 


The quantities with an asterisk are taken from the work: Nayler, Witing (1954), the remaining other 
experimental data from Kuhn, (1949). 


bonds of the molecule. Nevertheless, in the AS FEMO scheme for Anax 
get values closer to the experimental results than in the simple free-electron model. 


we generally 


This is mostly true for the singlet positions, which lie at about half way between the 
experimental value and that taken from the simple model. 

Kuhn (1949) has obtained good agreement with experiment for polyenes by 
introducing in the simple scheme a sinusoidal potential along the molecular core. 
In this picture the potential wells corresponded to z-electrons placed on the double 
bonds. The present results give evidence to the fact that the potential amplitude intro- 


TABLE II 


The theoretical locations of absorption maxima for symmetrical ‘cyanines 


crews 
ŽAS FEMO ÀpE ÅAS FEMO gg 
N sing. & tripl. simple sing. & tripl. simple 
arith. mean model arithm. mean model 


rj 
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duced by Kuhn should be smaller, і.е. the localization of electrons in not as strong as 
that previously considered?. 

Agreement between the theoretical calculations and the experimental data for 
the symmetrical cyanines constitutes another problem. This will be discussed in more 
detail in the following sections. It should now be emphasized that in comparison with 
the simple model the wavelengths of the absorption maxima are now shifted, as in 
the case of polyenes, in the direction of the shorter waves. The vinylene shift is less 
(approx. 1050 A instead of 1300 A) than for a simple model, its value, however, remains 
constant even for large N. Such behaviour is seen from the experimental data for the 
spectrum of symmetrical cyanines, and agreement between the theoretical and 
experimental numerical shift values is now аррагепз. 

The results of the application of the AS FEMO scheme may be refined by including 
the configurational interaction. It follows, however, from the calculations for ethylene 
and butadiene (see footnote 1) that the effect of this interaction is rather small. Never- 
theless, it should be expected that it may increase with the chain length, i.e. with an 
increase in the number of z-electrons. At the same time the terms occurring in the 
E(*3 ©.) — E(! Ф) which consist of Coulomb and exchange integrals will begin to 
play a more important role. 

It is characteristic for the calculations that the triplet levels fall almost exactly 
at the same place as the energy levels of a simple free-electron model. There seems 
to be no explanation for this fact at present. | 


$ 2. The Concept of the Resonance Barrier 


In the symmetrical cyanine cation two principal elements may be distinguished: 

1) —a chain of CH groups connected in turn by double and single bonds. The 
nitrogen atoms are placed at both of its ends. 

2) —а “nucleus” composed of the system of those atoms at one of the ends of 
the molecule, which do not form a chain of conjugated double bonds, and of those 
atoms of the chain which are directly bound with this system. As a rule from 1 the 
nitrogen atom and often one or more carbon atoms will belong to the nucleus. Symme- 
trical cyanine evidently has the same nuclei at both ends. 

There still exists a radical at the nitrogen atom. We exclude it from the system 
which we understand as the "nucleus". 


* However, the free-electron model can be applied with quite good results to polyenes with a small 
N/2 
number of z-electrons. Then the z-electron density in the ground state (22 v2 clearly corresponds to 


f n 
the concentration of these electrons in double bond regions (see Fig. a); for hexatriene see Platt (1954). 


With increasing chain length the maxima become flattened, we obtain a uniform density distribution, 
which corresponds to symmetrical cyanines. At the same.time it is interesting to note that the “experi- 
mental" potential amplitude given by Kuhn also decreases with the chain length. 
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The symmetrical cyanines are characterized by the equal probability that in the 
conjugated chain between any two atoms there exists a double or single bond. It follows 
that the charge is uniformly distributed along the atomic core, i.e. there is no collection 
of charge in the space corresponding to the bond of the pair of neighbouring atoms 
when the charge on the other bond is decreased. 

We would not, however, venture to say that the charge distribution is the same 
in all parts of the core. The end groups interact with the chain and the effect of this 
interaction will be, as we shall see, particularly strong in the regions of the conjugated 
system which are included in the nuclei. Outside the nuclei we shall assume, further, 
that the chain is approximately “unperturbed” and the charge is distributed uniformly. 

The electrons of the conjugated system in part of the cyanine nucleus are compo- 
nents of the electronic structure of the nucleus, as well the electronic structure of the 
chain. The nucleus forms a molecular system generally quite different and isolated 
from the rest of the chain. Thus, it often happens that the part of the chain correspond- 
ing to the nucleus plays a more important role in the structure of the nucleus as 
a whole than in the structure of the conjugated system. 

Now the electronic structure of the nucleus tends to take on the lowest energy 
state and the actual charge distribution will no longer be represented by the chemical 
formula of the nucleus, since this tendency of the electron assembly to form the most 
stable system has to be taken into account. In this way the nuclei act as if they forced 
the electrons of the conjugated chain to take part in their most convenient (i.e. from 
the energetical point of view) structure. 

Let us illustrate our discussion by the example: 


VN 


e С, en uec 
R R 
(I) 


Peet) 


| 
R Ч 


- In the former theory the actual charge distribution in the molecule is the inter- 
mediate value between (a) and (b). In the conjugated part of the nucleus this value 
lies between the state when both z-electrons are closely linked with the nitrogen atom 
and that, when they are distributed along the N — C bond forming a double bond. 
In the rest of the nucleus the charge distribution for both cases (a) and (6) remains 
unchanged; there are two conjugated double bonds. Together with a pair of electrons 


£ 
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from the conjugated part they form a system very similar to that of three double bonds 
in the hexagonal ring, i.e. the Kekulé structure of benzene. 

Such a system, as we know, is very unstable and reveals a strong tendency to admit 
lower energy corresponding to a uniform charge delocalization along the hexagonal 
core. Owing to their share in the electronic structure of the ring, the electrons from 
the chain tend to a state of much lower energy than that which would correspond to 
their participation in the conjugated system only. 


в 


The above considerations have their analogue in the theory of the resonance structures. In (I) two 
extreme resonance structures of a dye are shown. However, the actual state of the molecule cannot be 
obtained only from the resonance of the extreme structures. The interaction of the intermediate structures, 
corresponding to the position of positive charge at the C atoms in the chain, must also be taken into 
account. When the positive charge is placed on any of the nitrogen atoms we have a system of three 
double bonds conjugated in a ring. The stable benzene state is then more easily attained than in the case 
where there is no + charge on the М atom, i.e. when both electrons are in the closest neighbourhood 
of N. Therefore the limiting structures have lower energies than the intermediate structures. (See dis- 
cussion of Brooker (1942)). We shall further discuss this analogy later on. 


Let us now extend our treatment to the free-electron theory. It is well-known 
that the symmetrical cyanines seem to be best investigated by means of a free-electron 
model, because of the equivalence of the bonds in the chain.? 

In the box of the free-electron model we distinguish a part corresponding to the 
chain between the two nuclei and parts corresponding to the chain included in the 
nuclei themselves. According to the above discussion, in the first case we have, following 
Kuhn, the periodic potential, which for simplicity, we approximate by a constant value. 
Two deep potential wells caused by the attraction of the moving (л-) electrons by the 
end groups correspond to the other parts. Thus a potential barrier appears (see Fig. 1). 


60—0—0—€)—€)—€)—9) 


Fig. 1. The free-electron model with resonance barrier 


3 The decrease in the potential at the ends, which is caused by the presence of the nitrogen atoms 
alone, these atoms being more electronegative than the carbon atoms, is not considered here. This effect 
is rather small in comparison to the resonance barrier height, and does not select the cyanines. 
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Hereafter this will be called the resonance barrier, for, as we shall see later, there exists 
an analogy between our barrier and the intermediate structures in the resonance theory. 

We assume for simplicity that the barrier has a rectangular shape and that the 
ends of the barrier lie at the middle of the C — N bonds. This seems justified in most 
cases. The length of the free-electron path is equal to that in the problem without 
а barrier. 

Thus we have a quantization problem of free-electron motion in a box, where 
the potential at the ends is equal to zero, and a potential barrier of a constant height V is 
stretched over the carbon atoms. 

Assuming that at the particular sections of the box the solutions are given by the 
free-electron eigenfunctions* and taking the advantage of the condition of continuity 
of the functions and their derivatives, we get for the eigenenergies E, of the barrier 
problem the following equation?: 


Зл А ДЕ 2 а? \ 1% 
cos Tees sin 7 1 1 25 + 
. Зл 7 NO 2\% 2\ и 
еп b cos ee: 1 ( а) |. 1-5) cu (5) 


M TO = а? \ %2] a? 
+ sin? |—-———x]) sin | л: pale: 
+ sin (2 е sin Е Ар | P 0, 


Z (= N+ 1) is the number of atoms in the chain; 
h?a? ERE 
_8mL?’ 5. Вт 12 


It is now essential to find V. The value of V must be closely connected with the 


ТА 


(6) 


ability of the cyanine end groups to form stable systems. The measure of this ability 
is the energy difference of the unstable system, i.e. the nucleus without z-electrons 
of the chain as well as the -electron pair in C — N region, and the energy of a stable 
system, i.e. the nucleus as a whole. We shall call this energy difference the stabilization 
energy. It may also be obtained in a rough approximation on the basis of the free- 
-electron model. 


The energy of a z-electron pair in the C — М region is greater than that for a common double bond 
We assume that it corresponds to the energy of a pair of electrons in a box two C — C bonds long (instead 
of tree bonds for a common double bond system as in the case of ethylene). E 


For (1а) we divide the ring forming the nucleus into two segments: 1° — the 
free-electron path for C — N; 2? — the path in a system of two conjugated double 
bonds. Thus the stabilization energy is given by 

E DE Ey УР ERE (7) 


sta 


* The electronic interaction is neglected here. 

5 An equivalent eqvation was considered by Kuhn (1951) in his discussion of the efect of the 
electronegativity differences of the core atoms. However the present form of equation is better 
suited to our further discussion of the energy states. 


rj 
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For cyanine with a nucleus of the type 


: S 
"m t 3 3 $ ae (П) 
Bey een PS 123 
Y За pA d $e Q v DE 
N N 
| | 
R R 


the stabilization energy corresponds approximately to the difference of energies for 


a ring and a C — N segment 


BS 
D 
and the system 
Toe (IIb) 


where the free path is also lengthened by опе bond on the last atom, i. e. we get a system 
of styrene type. In addition, besides the (II)-scheme an alternative system is also 
possible: 


ii C 
T 3 a» a dl NA (IIc) 
УААН a 
| | 
R R 


therefore the stabilization energy for this nucleus ought to be calculated twice. 
The singular stabilization energy of the above type is realized approximately by: 


ya ju ys | 
ES 2 s pd 
A D pk P. Ха 

ПЕГЕ et М-Н О NERI (Ш) 
Sy n eer 


where it is assumed that the interaction of the conjugated chain with the ring through 
the C — С bond is weak and that only the z-electron pair from C — N region takes 
part in the stabilization of the nucleus. In fact, we see that for the same number of 
л-еЇесігопѕ (М = 10) in the conjugated chain there is a large difference in the wave- 
lengths of the absorption maxima in cyanines (II) and (III): for (II), where there 
is a higher barrier, this length is greater. _ 


On the Free-Electron Theory of Absorption Spectra 919 


When in cyanine nucleus (in the non-conjugated part) there are no additional 
double bonds (e.g. thiazoline system) then the end groups are weakly stabilized. The 
barrier is low and caused mainly by the presence of the N atoms in the chain. 

The stabilization energies for various cyanine nuclei are given in Table VIS. Later 
on they will be treated as V in the barrier problem. 


§ 3. The Effect of the Resonance Barrier on the Energy States 


We shall discuss the effect of the barrier of a variable height V (or а?) and of 
different lengths corresponding to Z = 5, 7, 9. From the numerical solutions of (5) 
it follows that x increases proportionally to V (see Fig. 2). This linear dependence 
is different fora particular л and Z, and is well satisfied for not too high V. For larger 
V the straight lines are slightly bent, forming approximately arcs of parabolas with a? 


е 
4 ee ажа) 


pt 
Beco, 


Fig. 2. The numerical solutions of Eq. (5) for the case Z = 5 (N = 6) and а?< 16. 
The origins of the straight lines correspond to the case when the barrier does not exist (a? — 0). We 
see directly from (5) that x is then equal to n from (1), i.e. we obtain the states of a simple free-electron 
model. The X are associated with the exact solutions 


as the axis. In addition it appears that at a fixed Z (a constant barrier length) the 

directional coefficients of the lines corresponding to increasing n's — decreases. Also 

the parabolas tend to meet. (For numerical results see Tables III and IV). 

Е ere wt era NE NITED how? nort pees et Tui Ins que got Teese ае Aer 
6 For the nuclei given in Table VI, which have not been considered above, we calculate the approxi- 

mate stabilization energies from the expressions: 


for number 4 — 2 [Ес N zs Eethyl. syst. Ebutad. cere 

for number 5 — (Estyrene syst. sr Ec. Nee Epapht. syst. ae Evenz. syst. F Ec -N7 Е; styr. syst. 

The values for the eigenenergies of some systems were taken from Scherr (1953). For energy calcu- 
lations the mean length of the bond was taken equal to 1.39 À (as for ring systems). 


£ 
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TABLE LE 
The — m — directional coefficients for x = та? + n lines corresponding to the approximate solution 
of eq. (5) for various п and Z 
The linearity is preserved well enough for Z = 5 up to а? = 16; for Z = 7 up to 2225: for —= 9 
up to a? = 36. Further on the lines for greater а? change into parabolas. 


N 
NS 2 3 4 | 5 6 
VA S | "a АР 
5 0.1205 | 0.0831 0.0683 | 0.0526 0.0383 
7 | 0.0917 0.0715 | 0.0620 0.0534 
9 | 0.109 0.0769 | 0.0631 0.0559 
TABLE IV 


The — b — and — c — coefficients in parabolic eq. а? = bx? + c (х > 0) corresponding to the approxi- 


mate solution of eq. (5) for greater а”. (See explanatory notes to Table III). 


== 


SET 
b pS 3 4 5 6 
A — —— 
5. | 1.217 | 1.767 
7 1.203 1.424 
9 1.142 1.240 
E ; 
c N 3 4 5 1 
Е UM | | 
Ет 045 — 31.12 
7.4 — 16.70 —81.21 
9 | | — 24.19 —43.15 


Thus, the lower the given quantum state, the greater, in general, the increase in 
energy for the given barrier height as compared to the corresponding eigenenergy 
of the nonbarrier problem. 

On the other hand it follows from the calculations that at à fixed height of the 
resonance barrier, and for a given n, x (and subsequently E.) increases with the chain 
length. This may be treated as a parallelism to Moffitt's conclusion from his analysis 
of absorption spectra of symmetrical cyanines on the. basis of resonance structures 
theory (1950). With an increase in the number of CH — CH groups in the chain, 
the number of intermediate structures ("bridge" structures) is increasing, and the 
differencé between the given state for the case when these structures are not considered 
and the state obtained by introduction of all the structures also increases. 

For a constant barrier height and for very large Z, eq. (5) gives solutions corres- 
ponding to results obtained by Kuhn for the constant potential over the entire box. 
of afree-electron model. This potential causes a shift of all the energy levels by a constant 
value V. Thus for very long chains the effect of the resonance barrier with a fixed height 
on the spectrum ceases to play essential role. 
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In general, if the height of a resonance barrier is known, the wavelength of the 
appropriate absorption band may be calculated (Table VI), and vice versa. 

We shall evaluate for the latter case the height of the resonance barrier for various 
chain lengths starting with the experimental data of the spectrum of cyanine: 


a 


; S 
aa aS А zt p e 
| R С (CH = СН), сн с d 
х 4 N A и 
мы. NÁ V (IV) 
| | 
R R 


Thus we obtain: 


І 2 (А) barrier A (А) 
exp. height (eV) | simpl. model 

0 4230 | 11.6 | 3320 

1 5515 | 11.3 | 4595 

2 6500 | 10.1 | 5880 


and the theoretical barrier height (the stabilization energy) lies above 9.8 eV. Hence 
it follows that with an increase in chain length, the barrier corresponding to the parti- 
cular (ТУ) homologues subsequently decreases. This may be explained perhaps by 
the growth of the self-stabilization of a cyanine chain with an increase in the length 
of this chain. The end-group effect (as a relative stabilization effect) is smaller for the 
longer systems than for the shorter ones. à 


We may regard this result as being experimentally proved by the fact that the ability to bind a proton 
to the cyanine (IV) chain increases with an increase in the number of CH — CH groups (Brooker 1942). 

When the proton is joined to the chain, it spoils the delocalization of z-electrons and we get a system 
partly similar to the model for polyenes. As mentioned before a periodic potential corresponds to such 
a model (Kuhn). Then the mean value of the potential is close to zero, so the potential barrier is destroyed 
here. We could expect that the facility with which the barrier is destroyed is connected with the barrier 
height. 


§ 4. Conclusions 


When the free-electron model is adapted to the treatment of absorption spectra 
of symmetrical cyanines, a question may arise as to what extent this model fulfils its 
purpose. 

The old simple FE model, developed for cyanines by Kuhn, has given good agree- 
ment with experiment only for a few molecules of this kind. For the majority of cyanines 
the A,,, experimental values were less or greater than expected by the model; for 
the latter case the differences used to be very large. The model did not explain the 
essential question: Why the conjugated systems in cyanines of the same number 
of -electrons often gave big differences in the absorption wave-lengths? 
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Kuhn (1948, 1949) tried to explain the absorption in the longer wave-length 
region by assuming an additional lengthening of the free-electron path, but he did 
not give any theoretical estimate of it. For the shorter wavelengths one should have 
then assumed a shorter path than that for a free-electron in the molecule. Now, by 
а more precise treatment of the AS FEMO method all the experimental values for 
Anax lie in the region of wavelengths longer than those predicted from the model. 
However, in the case of molecules for which the physical conditions of a pure free- 
-electron model are best satisfied, A,,,, should be in quite good agreement with the 
AS FEMO values. 


In fact, for molecules (Simpson cyanine dyes) 


ax 


MeN = CH — (CH = CH), — N Me? (V) 


in which there is no stabilization in the end-groups, we obtain: 


TABLE V 
The Simpson cyanine-dyes 
IM a А 
j А AS FEMO Experim. 
E 

0 6 2900 3090 

1 8 3905 4090 

2 10 4980 5110 


When we introduce the resonance barrier into the simple free-electron model, 
we can roughly, at least, overcome the difficulties which arise in the interpretation 
of the spectra by а non-barrier model. As follows from the solutions in Sec. 3, the 
presence of a high barrier raises the ground state more than an excited state. Thus 
the difference between these states decreases and we get absorption at longer wave 
lengths than for the non-barrier model, as is observed. This change of the position 
of states is determined, as we have seen before, by the barrier height, and therefore 
it is a function of the stabilization energy of the cyanine nuclei. 

In our picture the evident shift towards the long waves appears, as a rule, for 
greater barrier heights. We very often deal with cyanines having such high barriers. 
For lower barriers the results are generally not much different from those obtained 
on the basis of a simple model. We even get absorption in wavelengths shorter than 


* In addition to Kuhn, the following authors, on the basis of empirical data, also took into account 
the effect of the end nuclei: Nikitine (1951), Basu (19542) fitted the length of the free-electron path; Fukui 
et al. (1953) fitted the periodic potential in the box, Simpson (1948) fitted both the C — C bond distance 
and the path length. The non-empirical treatment of Basu (1954b) is rather doubtful in view of the assy- 


metry of the potential in the box and very good results seem to be obtained for an exceptional case 
only; See, however, footnote 9. ` 
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predicted by the model, the deviations, however, are small. Thus we may conclude 
that, in general, for barriers that are not high the /,,,, are close to the А from the non- 
-barrier model, as is actually the сазе8. (Table VI). 


Fig. 3. The shift in the energy levels of the ground state (g.l.) and the excited state (e.l.) caused by 
the existence of a high resonance barrier. 
(a) free-electron model without-a barrier 
(b) a barrier model 

The change in the energy corresponding to the absorption (A £) is significantly greater in (a) than in (5) 


We shall now discuss briefly the vinylene shift problem. Let us calculate the 
wavelengths of absorption maxima for a fixed barrier height in the case of several 
subsequent homologues produced by the increase in the number of vinylene groups 
of the chain. We find that the change in À, 
homologues (the vinylene shift), particularly at the high barriers, is much greater when, 
for instance № = 6 -> № = 8, than for №=8 > № = 10. For the transitions corres- 
ponding to still greater N this shift approaches that given by the non-barrier free- 


max corresponding to the two subsequent 


-electron model. 

For ^?" cyanine, for example, the first shift calculated from the barrier model is greater by 150 À than 
the next one. (Table VI). The simple model does not show any difference in the shift here. 

Thus we obtain the theoretical explanation of the well-known experimental fact 
that the first vinylene shift for the given series of cyanine is greater than the subsequent 
values of this shift.? 


8 From the numerical point of view the differences between the results of the barrier model and 
the experimental data sometimes still remain serious. On the other hand, the roughness of our method is 


apparent. To improve it the effect of the radicals on the ends and the difference in the electronegativities 


between C and N should be taken into consideration in addition to the theoretical refinement of the 
calculations. 

9 Such behaviour of the vinylene shift was first communicated by Brooker (1942). He listed cyani- 
ne “3“ as a typical compound (Table VI). This result is not confirmed, however, by the experimental 
data of cyanine “5“. Here the first shift is smaler than the two next, which are approximately the same. 


_ (These data are not yet very certain in view of the fact that location of the absorption maxima is not 
_ clear — see Мите (1951)). On the other hand the height of the resonance barrier suggested by the 


experimental data (as well as by calculations) is extremely large (the greatest shift of absorption maxima 


- toward long waves). Then the probability of leakage" of electrons across the barriers at the ends of the 


| 


box outside the nuclei approches that across the resonance barrier. This suggests another treatment 


of the problem. M 
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TABLE VI 


'The UTES energies and location of absorption м о barrier model 


Stabilization 1 (А) A (À) 
The Е 
energy (eV) teor,. exp. 


nu йг: к: SIRE Я 
| 2 v6, Bis E 
s ets " > 9 
- ne L 
Now eaa vier И 
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TABLE VI (continued) 


Stabilization 


Nr The nucleus 
energy (V). 
SS 
1“ | | 
а, И 
Кт N 6 11.8 4300 5235 
R 
8 11.78 5725 6040 
10 11.78 6925 7080 
z 12 11.78 8265 8200 


Only one of the extreme structures of the nuclei is shown in the table. For the experimental 
data see Brooker (1942), Brober etc. (1941), Kuhn (1948) Nikitine (1951). 

For large NV the theoretical results of the barrier model are quite close to those of the 
non-barrier model. It is easily seen, however, that the barrier model as a whole approaches con- 
siderably better the experimental data than the simple one (see Table IT). 


In Sec. 2 we indicated the analogy which existed between the barrier in a free- 
-electron model and the resonance formalism. Let us now compare with the results 
of this work Moffitt's (1950) conlusions from his detailed analysis of electronic states 
of the symmetrical cyanines on the basis of the resonance theory. We see that in many 
points the effect of intermediate structures on the spectrum is similar to the effect 
which is caused by a barrier in the free-electron scheme. 


'Thus we have: 


1) Moffitt: "the excited structures have greater effect on the ground state than 
on the excited states". 

The present theory: "the introduction of the barrier causes a greater change in 
a ground state than in an excited level. The higher is the level, the smaller is the effect 
of the barrier". 

2) Moffitt: "with an increase in the chain length the energy difference of the state 
corresponding to the extreme structures and the state obtained from the interaction 
of all the structures increases". 

The present theory: "with an increase in the chain length at a fixed barrier height, 
the energy difference of the state of a non-barrier problem and the same state in the 


presence of the barrier increases". 
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3) Moffitt: "with the increase of the chain length the mutual stabilization of the 
ground state and the excited state decreases. The vinylene shift is, therefore, non- 
-convergent”’. 

The present theory: "with the increase in the chain length (and consequently 
with the increase in n) at a fixed barrier height the situation becomes similar to that 
in the non-barrier free-electron model. This is known to causea constant vinylene shift." 


КРАТКОЕ СОДЕРЖАНИЕ 


C. Ольшевский Теория абсорпционных спектров некоторых ‘линейно сопряжённых 
систем 


В настоящей работе применён метод антисимметрических молекулярных 
орбит свободного электрона к вычислению абсопционных максимум в спектре 
симметрических цианинав и полиенов. Результаты продискутированы в сравнении 
с эмпирическими данными и результатами простой модели свободного электрона. 
Потом показано, что можно получить значительное улучшение теоретических 
ценностей для цианинов посредством введения в модель свободного электрона 
резонанского барьера потенциала. Высота барьера определена энергией стабили- 
зации ядер цианина. Подробно продискутировано влияние на энергетическое 
состояние молекулы, вызнавнное введением в простую модель свободного элек- 
трона потенциального барьера. Указано на согласованность полученных зтим 
путем результатов с выводами теории резонанса. 
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DOMAIN FORMATION IN FERROELECTRICS 


By ARKADIUSZ JAŚKIEWICZ 


Physical Laboratory B. Bierut University, Wrocław 


(Received September 22, 1956) 


The dielectric properties of ferroelectrics depend on the domain structure and 
domain dynamics. The formation and movement af domains have been studied by 
various authors: Forsbergh (1949), Merz (1954), Little (1955). Our main aim is to find 
how the formation of the c domain of BaTiO, depends on the electric field. 

The nucleation theory supposes that local fluctuations in free energy produce 
small, ordered regions called embryos. The free energy has to overcome various energy 
expenditures. Owing to the increase in the surface energy and other expenditures 
associated with the formation of an embryo, most embryos are unstable and vanish. 
Embryos greater than a critical size are stable and grow with a decrease in the free 
energy of the substance. This implies that a critical size can be reached only by the 
kinetic process of statistical fluctuations see e.g. R. Smoluchowski 1951. 

The change in free energy AF associated with the formation of a nucleus in ferro- 
electrics can be represented by the expression 


AF = ДЕ, + AF, + AF, + AF, 7 (1) 


The first term represents the surface energy (wall energy); the second, the energy 
volume energy (dependent on the parameter of order); the third, the energy of inter- 
action between the polarized dielectric and the field; the last term, the depolarization 


energy. 
If we take the spontaneous polarization as the parameter of order, then for the 


order energy we obtain the expression 
AE, = аР? (2) 
where æ depends on the shape of the embryo, V is the volume, and Р, is the spontaneous 


polarization. 
| | , (227) 
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Formüla (1) giving the free energy change in the absence of a field, can be written 
in the form: 


AF = oS aP V + ВРАГ (3) 


where д is the wall energy/cm?, S the wall area, В the depolarizing factor. 

If a phase transformation is to occur, x must be greater than f, because the free 
energy of the new phase is lower than that of the metastable phase. 

We assume that the embryo has a cone-like shape with length / and radius r in 
accordance with optical observations Hippel (1950), Merz (1954). Using (3) we have 


1 
ДЕ = ло (r? + rl) (a p) 3 лг?ІР? (4) 
if we write 
a= p =y (5) 
we get 
AF = no(r? + rl) — y = nr?lP? (6) 


We obtain the critical dimensions from the maximum of (6): 


CRY 7 Oo; z 
e ie R 
We see that there will be only large domains about the transition point. The 
nucleation rate of domains we can get from well known expression for the rate of 
nucleation. 
Now, we take under consideration a ferroelectric crystal in an external field PF. 
The free energy change of nucleus formation is given by 1, thus, from (6) we can write: 


AFP = xo(r? + rl) — = ynr*lP? — EP,V (8a) 


AF* = mo(r? + rl) — E yartlP? + EPV (8b)- 


where — EP,V is the electrostatic energy of a nucleus with parallel polarization, ЕР, V 
the energy of a nucleus with antiparallel polarization. 

If we repeat our procedure we obtain the critical dimensions: 

a) domains with parallel polarization 


3c 6c 
ГР а ; P = : 
РОР, Е)‘ * ~ POP, +5) us 
b) domains with antiparallel polarization 


га = See Tee a e. n e 6c 
PPLE) infissi tpa (2b) 
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We see that equations (9a) and (9b) explain the behaviour of domains in the electric 
field. In the case of a nucleus with antiparallel polarization, the critical dimensions 
are greater than in the case of parallel polarization. All the antiparallel nuclei which 
have dimensions less than the critical one are unstable and vanish. The destroyed anti- 
parallel domains have a greater probability to nucleate parallelly than antiparallely. 
The probability may be obtain from the nucleation rate. The last equations are well 
suited for a study of domain formation in connection with experiments such as those 
of Piekara and Pająk 1953 and of Badian and Jaśkiewicz. 


The author wishes to express his thanks to Dr. Z. Galasiewicz and Professor B. Ma- 
kiej for helpful discussions. 
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NOTE ON THE ??C(d,p)'?C REACTION NEAR THE 4 MeV RESONANCE 


J. DaBRowsk1 AND B. TULCZYJEW 
Institute of Theoretical Physics, University of Warsaw, Warsaw 
(Received February 11, 1957) 


Recently several observations have been performed for the stripping reactions 
in the vicinity of an isolated resonance of the compound nucleus. A calculation scheme 
for these reactions is outlined in the present note in its most simple version for the case 
of the 12С(4,р)13С reaction near the 4 MeV resonance investigated experimentally 
by Bonner et al. (1956a).! We adopt their values: J = 3* and the width Г, = 0.15 MeV 
of the excited level in !4N. 

The formula for the cross-section for the reactions of the type considered here 
has been given by Thomas (1955) (cf. also Dabrowski 1956?). To perform the numerical 
calculations on the basis of this formula we make the following simplifying assumptions: 

(1) The Coulomb effect is not taken into account. 

(2) The potential scattering of deuteron and proton is neglected. 

(3) Following the suggestion of Bonner at al. (1956a) we assume the values 1, = 
jp = 0 (channel spin) as responsible for the compound nucleus mechanism. Thus 
we neglect the contribution of 1; = 4 and j, = 1. 

(4) The neutron is captured with /, = 1, this being in agreement with the shell 
model. 

Owing to these assumptions we may write the differential cross-section for the 
1?C(d,p) 13С reaction in the vicinity of the 4 MeV resonance in the form 


do/dQ = (уа) [(do/dQ), + A(do/dQ), + 249149) с], 


where y? is the reduced width of the neutron bound state іп 13С; А = Veal V3 Vo Va 
are the corresponding partial widths of the virtual level in UN; y?a-l(da|dQ)s is the 
usual stripping cross-section; y2a~1A(do/dQ), is the interference cross-section (con- 
structed as in Eq. (17) of the paper by Thomas (1955); y?a~14?(do/d2) с is the Breit- 
Seren tee STR ee etal eee а: бое сы ЗО Ос ЕЕ ЕС ОК 
1 [t is to be noticed that in the paper by Bonner et al. (1956a) the energies given in Fig. 3 should 
be reduced by about 0.08 MeV and that the solid curve in Fig. 11 is not correct (it should have no minima 
for 0° and 180°). 
2 The symbols used here are the same as in that paper. 


, (231) 
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Wigner — one level — cross-section (without the Coulomb and potential scattering 
phase shifts). 

The calculations were performed with the help of Butler’s theory and Born’s 
theory (transparent nucleus). It was found that the best value of the (stripping) radius 
of the 12С nucleus is а ~ 6 x 10-13 cm (for the Born theory still a larger value of 
a should be used). 

Now let us discuss the value of the resonance energy. Owing to the observation 
of Bonner et al. (1956a) the deuteron resonance energy E% = 4.004 MeV, and 
the occurence of the resonance is most distinct in the first stripping maximum. The 
interference cross-section causes that the true resonance energy E, differs from E%. For 
small Я we get the maximum value of do|dQ in the first stripping maximum for 
E95 ~ Е, + Г/2. In our calculation the plus sign is used (E, = 3.93 MeV) since 
in the reactions 19B(a, п)13М and !B(a, p y)!*C (Bonner et al. 1956b) the corresponding 
resonance energy was observed to be equal to 3.98 MeV, and it is hoped that in the 
case of both these reactions the role of the direct interaction is smaller and the observed 
resonance energy is closer to the true resonance energy than in the case of the (d,p) 


reaction. 
The quantities A, y were determined by fitting the values of the theoretical cross- 
-section at the first stripping maximum to the experimental values at E; = 3.84 and 


4.01 MeV. We found 


0.87 x 10-13 MeV cm, 0.0099 MeV, for the Butler theory, 
== оа m 

0.22 х 10-13 MeV cm, 0.0012 MeV, for the Born theory. 
The corresponding contribution of (0/40) с to do/dQ at the first maximum is about 3% 
for both the energies. It should be noticed that the dimensionless reduced width 
0? = 2M*ay?/h? in the case of the Born calculation (0.058) is in rather close agreement 
with the value 0.07 given by Thomas (1952). 

The corresponding angular distributions are given in Fig. 1. It is hardly possible 
to decide whether the Butler or the Born angular distribution is better. We see from 
Fig. 1. that the compound nucleus mechanism is not strong enough to give the observed 
hack scattering maximum. This is still more evident from Fig. 2. for E, = 4.12 MeV. 
On the other hand, for this energy the compound nucleus mechanism is too large, 
since the calculated value of the cross-section at the first maximum is larger than the 
experimental value. 

Assumptions (1), (2) are very crude. One may argue that the influence of the 
Coulomb effect and potential seattering has been investigated by several authors 
(Horowitz and Messiah 1953, Yoccoz 1954, Grant 1945—1955, Tobocman 1954, 
Tobocman and Kalos 1955). But the neglect of these effects in the interference cross- 
-section is a far reaching simplification, since it follows from it that (do/dQ ) ~ (E,—E,). 
Of course, the term proportional to (E; — Е) will remain in a more exact formula 
also and will cause a shift of the observed resonance. 
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Assumption (3) is somewhat arbitrary. It is to be noticed, however, that ап intro- 
duction of other angular momenta into the calculation would introduce new parameters 
(the corresponding partial widths) automatically?. This would destroy the simplicity 
of the present treatment. 

The results of Bonner et al. (1956a) suggest the vanishing of (do/dQ), near the 
first stripping maximum (cf. their Fig. 7). This is not the case in our calculation. Con- 
nected with this is our rather unrealistic assumption E, = 3.93 MeV. All these facts 
show that our simplifying assumptions are going too far. The spin value 3+ of the re- 
sonance level is rather arbitrary also. Nevertheless it is hoped that one may conclude 
already on the basis of this calculation that the compound nucleus formation is not 
a sufficient reason for the back scattering maximum even in the vicinity of the reso- 
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Fig. 1. Angular distribution of protons from 1C(d,p)!°C for the deuteron energy E; = 3.84 and 4.01 MeV 


3 The situation would be simpler in the case of spinless particles, where the number of the possible 
combinations of the angular momenta becomes smaller. For such cases the scheme presented here would 


be more useful. 
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nance. Further it is hoped that the scheme presented may be of some value for other 


similar calculations. | 
We are indebted to Dr. J. В. Marion for explanatory remarks concerning his 


papers (with Bonner et al.) and the problem discussed in this note. 
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Fig. 2. Angular distribution of protons from !?C(d,p)!3C for the deuteron energy Ej; = 4.12 MeV 
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Derivation of the Final State of the Outgoing Neutron Eq. (3) represents the 
"equivalent" outgoing state of the neutron scattering in a spin-orbit potential its 
radial partial waves being proportional to J, (k, г) — В (Г, Л) АФ (k, r). Eq. (3) 
results as a generalization of Toboeman’s Eq. (2) for spinless particles (1954, Eq. 
(19)). In analogy to Tobocman's Eq. (2) for spinless particles let us introduce the 
spin-dependent function: 


ey A) ТОТО (1’) 


W being the exact three-body function of the problem. 
In analogy to Tobocman's Eq. (11) we may write the integral EI fort (15: 


Е М itn | "н" | "Een Р 
о dr. Sgr ae Yi V "n y CE On) (2) 
Sug (n> Or) 2лћ? i м es 


У, being an integral operator. We can expand fur (ras On) in total angular 


momentum eigenfunctions: 
Ж л My 5 
Fug (r,, с.) mE >) JH (r,) UE ji (Q,, On) (3°) 
LM, УМ; 3 


where: 


a 

C (Qas о) = У. (25 ИД му LMj, (Q,) X,; (o,) (4) 
Mp, Hy 

Introducing the expansion (3') into Eq. (2) gives as a final result: 


2k. M xy M; , 7, 
"cs (r, ) = ih? dr, do, Л. (А, г2) hy (5, Г) Ч (©, 2 Op) V pu ur (r, ? on) 


(5) 


(235) 
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where: 
"x5 TOTEM ip, NATI ТОСТИ ars 
yi. fom "nga for r, <r 
la or Tn n Tn n п 


One can now еаѕШу eliminate the interior region of the nucleus and to derive 
the expression for the distortion coefficient f//(L, J) in complete analogy to 
Tobocman's expression for fj; , м (Eq. (17)): 


д д 
дель, — Л) 5; Im ЛЮ 


9, 
дг 


p" (Г, J) = | (6) 


MP (6,7) — MP (5,7) © In ff () 


If now, we indentify F In р o with the logarithmic derivative appro- 
Д R, 
priate for free neutrons incident on the residual nucleus, the и,-дерепдепсе is 


neglected and we obtain В (Г, J) = 22200 > пу being the average reflection 


coefficients in analogy to Eq. (IL5) — (П-7). Eb 
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